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ABSTRACT 



> 

(^-^ Aims. By using the full span of multicolor data on a representative sample of first overtone RR Lyrae stars in the Large Magellanic Cloud 
QQ (LMC) we revisit the problem of the incidence rate of the amplitude/phase-modulated (Blazhko) stars. 

■ Methods. Multicolor data, obtained by the MAssive Compact Halo Objects (MACHO) project, are utilized through a periodogram averaging 
Cn| ' method. 

Results. The method of analysis enabled us to increase the number of detected multiperiodic variables by 18% relative to the number obtained 
by the analysis of the best single color data. We also test the maximum modulation period detectable in the present dataset. We find that 
' variables showing amplitude/phase modulations with periods close to the total time span can still be clearly separated from the class of 
, stars showing period changes. This larger limit on the modulation period, the more efficient data analysis and the longer time span lead to a 
^ I substantial increase in the incidence rate of the Blazhko stars in comparison with earlier results. We find altogether 99 first overtone Blazhko 
stars in the full sample of 1332 stars, implying an incidence rate of 7.5%. Although this rate is nearly twice of the one derived earlier, it is 
still significantly lower than that of the fundamental mode stars in the LMC. The by-products of the analysis (e.g., star-by-star comments, 
distribution functions of various quantities) are also presented. 
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1 . Introduction stars. In a similar study on a very large sample containing 71 10 

RR Lyrae stars in the LMC, Soszynski et al. (2003) derived 
Current analyses of the databases of the microlensing 15% and 6% for the RRO and RRl stars, respectively. In ear- 
projects MACHO and OGLE (Optical Gravitational Lensing j^^^ ^^^^-^^ ^„ ^j^^ ^^^^ g^l^^y^ ^^^^^ ^j^^ observations of 
Experiment) have led to a substantial increase in our knowledge ^j^^ MACHO project, Alcock et al. (2000, 2003, hereafter AOO 
on the frequency of the amplitude/phase-modulated RR Lyrae ^^3^ respectively) got rates of 12% and 4% for the above 
(Blazhko [BL] ') stars. Moskalik & Poretti (2003) analyzed ^j^^^^^ variables. One may attempt to relate these inci- 
the OGLE-I data on 215 RR Lyrae stars in the Galactic Bulge. ^^^^^ ^^^^^ ^j^^ metallicities of the various populations (e.g.. 
They found incidence rates for the fundamental (RRO) and first Moskalik & Poretti 2003), but the relation (if it exists) is cer- 
overtone (RRl) Blazhko stars of 23% and 5%, respectively. A ^^^^^ ^ ^-^^^^ (^^^.^^ 2005; Smolec 2005). 
similar analysis by Mizerski (2003) on a much larger dataset 

from the OGLE-II database yielded 20% and 7%, respectively. Except for the SMC, all investigations indicate much lower 
These rates, at least for the RRO stars, are very similar to the incidence rates for RRl-BL stars than for RRO-BL ones. It is 
ones suggested some time ago by Szeidl (1988) from the var- believed that the cause of this difference is internal, i.e., due to 
ious past analyses of rather limited data on Galactic field and real difference in physics and can not be fully accounted for by 
globular cluster stars. One gets smaller rates by analyzing the the possibly smaller modulation amplitudes of the RRl stars. 
RR Lyrae stars in the Magellanic Clouds. Based on the OGLE- Since this observation may have important consequences on 
II observations of a sample of 5 14 stars in the Small Magellanic any future modeling of the BL phenomenon, we decided to re- 
Cloud (SMC), in a preliminary study, Soszynski et al. (2002) analyze the MACHO database and utilize all available obser- 
obtained the same rate of 10% both for the RRO and RRl vations (i.e., full time span two color data). In Sect. 2 we sum- 
marize the basic parameters of the datasets and some details 



' See Szeidl & Kollath (2000) for a historically more precise possi- of the analysis. Section 3 describes our method for frequency 
ble nomenclature. spectrum averaging. The important question of the longest de- 
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Table 1. Properties of the LMC RRl datasets analyzed in this 
paper 



Set 


(Tm) 




Colors 


#1 


6.5 


700 


'r', 'b' 


#2 


7.5 


900 


'r', 'b' 



Notes: (rtot) = average total time span [yr]; (A^d) = average number 
of datapoints per variable; Colors: MACHO instrumental magnitudes, 
see Alcock et al. (1999) 



tectable BL period from the present dataset and the concomi- 
tant problem of variable classification are dealt with in Sect. 4. 
Analysis of the RRl stars with their resulting classifications 
will be presented in Sect. 5. Finally, in Sect. 6 we summarize 
our main results with a brief discussion of the current state of 
tiie field. 



2. Data, method of analysis 

For comparative purpose, in the first part of our analysis, we use 
the same dataset as the one employed by AOO. Our final results 
are based on the full dataset spanning ~ 7.5 years. Basic prop- 
erties of these two sets are hsted in Table 1 . Both sets contain 
the same 1354 stars and cover the fields # 2, 3, 5, 6, 9, 10, 11, 
12, 13, 14, 15, 18, 19,47,80,81 and 82, sampling basically the 
LMC bar region. Because the selection of the stars was made 
earher on simple preUminary criteria such as period and color, 
some variables, other than first overtone RR Lyrae stars, were 
also included. Fortunately, the number of these other variables 
is only 22, that is small relative to the full sample. 

We employed a standard Discrete Fourier Transform 
method by following the implementation of Kurtz (1985). All 
analyses were performed in the [0, 6J d"' band with 150000 
frequency steps, ensuring an ample sampling of the spectrum 
line profiles even for the longest time series. The search for 
the secondary frequencies was conducted through successive 
prewhitenings in the time domain. The first step of it consisted 
of the subtraction of the main pulsation component together 
with its harmonics up to order three. In all cases the frequency 
of the given component was made more accurate by direct least 
squares minimization. When the two colors were used simul- 
taneously, we allowed different amplitudes and phases for the 
two colors and minimized the harmonic mean of the standard 
deviations of the respective residuals. Although somewhat ar- 
bitrarily, we considered the harmonic mean as a useful function 
for taking into account differences in the data quahty. 

For the characterization of the signal-to-noise ratio (SNR) 
of the frequency spectra we use the following expression 



SNR = 



(Ay) 



(1) 



where Ap is the amplitude at the highest peak in the spectrum, 
(Ay) is the average of the spectrum and cta, is its standard de- 
viation, computed by an iterative 4cr clipping. 



3. Spectrum Averaging Method (SAM) 

Unhke OGLE, where measurements are taken mostly in the 
Ic band (e.g., Udalski, Szymahski & Kubiak, M. 1997), the 
MACHO database has roughly equal number of measurements 
in two different colors for the overwhelming majority of the ob- 
jects (e.g., Alcock et al. 1999). This property of the MACHO 
database enables us to devise a method that uses both colors for 
increasing the signal detection probability. 

It is clear that simple averaging of the time series of the dif- 
ferent color bands cannot work, because: (a) the two time series 
are systematically different due to the difference in colors; (b) 
sampling may be different for hardware, weather or other rea- 
sons. Therefore, we are resorted to some spectrum averaging 
method that uses the frequency spectra of the two time series 
rather than their directly measured values. Obviously, SAM is 
also affected by the time series properties mentioned above. 
However, this effect is a weaker, because: (i) we are interested 
in the positions of the peaks in the spectra and these are the 
same whatever colors are used; (ii) spectral windows are sim- 
ilar, unless the samplings of the two time series are drastically 
different. As an example for the fulfilment of this latter con- 
dition for the MACHO data, in Fig. 1 we show the spectral 
windows in the two colors in one representative case. 

6.6094.5606 



•a 

3 



I 

•e 



O 



-1 1 

Freqtiency [d '] 

Fig. 1. Example of the similarity of the spectral windows of the 
MACHO instrumental 'r' and 'b' time series. Close-ups of the 
main peaks are displayed in the insets. The MACHO identifi- 
cation number of the star analyzed is shown in the upper left 
corner of the figure. 

In order to optimize noise suppressing, we compute the 
summed spectrum by weighting the individual spectra by the 
inverse of their variances 




cri + o-fy crt crt I 



(2) 



where Ri and B, are the amplitude spectra of the 'r' and 'b' 
data, (Tr and <Tb are the standard deviations of the spectra. 
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Before we discuss the signal detection capability of SAM, 
it is necessary to give significance levels for periodic signal de- 
tection when using the present dataset. Perhaps the simplest 
way of doing this is to perform a large number of numeri- 
cal simulations and derive an empirical distribution function 
of SNR. For this goal we generated pure Gaussian noise on the 
observed time base of 10 randomly selected stars. For each ob- 
ject we generated 1000 different realizations and computed the 
frequency spectra in our standard frequency band of [0, 6] d"' . 
The empirical distribution function was derived from the 10000 
SNR values computed on these spectra. Figure 2 shows the re- 
sulting functions for the single color and SAM spectra. For fur- 
ther reference, the 1% significance levels are at 6.9 and 6.5 for 
the single color and SAM spectra, respectively. All signals that 
have lower SNR values than these, are considered to be not de- 
tected in the given time series. The lower cutoff for the SAM 
spectra is the result of the averaging of independent spectra 
which also results in a greater number of independent points 
in the SAM spectra and concomitantly, in a slightly different 
shape of the distribution functions (see also Kovacs, Zucker & 
Mazeh 2002). 



A 

I 

CO 




Fig. 2. ProbabiUty distribution functions of SNRs of the ampli- 
tude spectra of pure Gaussian noise generated on the time base 
of the 'r' data (dotted line), 'b' data (dashed line) and their 
combination by SAM (soUd Une). 

In turning to the signal detection efficiency of SAM, we 
note the following. At very low SNR we expect no improve- 
ment by applying SAM, because the optimally achievable noise 
suppression is insufficient for the signal to emerge from the 
noise. In the other extreme, when the the signal is strong, the 
improvement is expected to be minimal and secure detection is 
possible without employing SAM. In order to verify this sce- 
nario, we conducted tests with artificial data. The result of one 
of these tests is shown in Fig. 3. The test data were generated 
on the 7.5 year time base given by the variable 10.3434.936. 
The signal consisted of a simple sinusoidal with a period of 
0.278 d and ampUtude A, where A was chosen to be the same 
in both colors and changed in 100 steps between 0.01 and 0.3 to 
scan different SNR values. At each value of the amplitude we 
added Gaussian noise of cr = 0.11 to the signal. This procedure 



was repeated for 100 different realizations. At each color and 
amplitude we computed the average and standard deviation of 
the SNR values obtained from these 100 reahzations. In order 
to characterize the gain by employing SAM, we also computed 
the ratio of the SNR values (^snr) as derived from the single 
color data and by the application of SAM. 
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Fig. 3. Uppermost panel: variation of the ensemble average of 
^SNR = SNRsAM/SNRcoior as a function of the test amplitude. 
Open circles refer to the 'r' data, crosses to the 'b' data. Middle 
panel: close-up of the low/mid SNR region of the above di- 
agram. Error bars show the Icr ranges of the values obtained 
by the various noise realizations. Lowermost panel: ensem- 
ble average of SNR. The dots are for the SAM values. Short- 
and long-dashed Unes show the 1% noise probability levels at 
SNR = 6.9 and SNR = 6.5, corresponding to the single-color 
and SAM spectra, respectively (see also text and Fig. 2 for fur- 
ther details). 



We see from Fig. 3 that there is a region where the signal 
has too low amplitude to be detected in the single color data, 
but not low enough to remain hidden in the SAM spectra. This 
is the most interesting amplitude/noise regime, where the ap- 
plication of SAM results in new discoveries. At higher ampli- 
tudes SAM leads 'only' to an increase in SNR. It is seen that 
the maximum increase in SNR is at around V2, as expected 
from elementary considerations. 

In order to show non-test examples for the signal detection 
efficiency of SAM, in Figs. 4 and 5 we exhibit the method at 
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Table 2. Number of single- and multi-periodic variables in set 
#1 of Table 1. 



Color 


r 


b 


SAM 




978 


895 


813 




376 


459 


541 


N2/{Ni+N2) 


27.8% 


33.9% 


40.0% 



Notes: N[ = number of single-periodic variables; N2 = number of 
multi-periodic variables 



work for the observed data of two variables. In both cases the 
single color data are insufficient for detection, but the combined 
spectra show clearly the presence of the signal. 



13.5842.2468 

SNR = 5.00 Ri 

1 - 




k , 

1 2 3 4 5 6 

Frequency [d ] 



Fig. 4. Example of the signal detection efficiency of SAM at 
low SNR level, when the single color data do not show the 
presence of the signal. Upper and middle panels show the fre- 
quency spectra for the 'r' and 'b' data, whereas the lower one 
exhibits the SAM result. The MACHO star identification num- 
ber is shown at the top, whereas the SNR values are given in 
the corresponding panels. All spectra are displayed in the same 
(arbitrary) scale. 



By analyzing all the 1354 variables of set #1 we can count 
the number of cases when significant components are found af- 
ter the first prewhitening. The result of this exercise is shown 
in Table 2. We see that there is a strong support of the obvi- 
ous assumption that application of two color data significantly 
increases the detection probability of faint signals. The gain is 
18% relative to the best single color detection rate. 



3.6597.756 

SNR = 5.50 Ri 



1 - 




k 

1 2 3 4 5 6 

Frequency [d ] 



Fig. 5. As in Fig. 4, but for the case when the signal is at the 
verge of detection in the single color data and SAM increases 
the rehabiUty of the detection. 

4. The phenomenological definition of the BL stars 

- Separation of the BL and PC variables 

The current loose definition of the BL stars relies purely on the 
type of the frequency spectra of their light or radial velocity 
variations. Usually those RR Lyrae stars are called BL variables 
that exhibit closely spaced peaks in their frequency spectra. 
Although most of the variables classified as type BL contain 
only one or two symmetrically spaced peaks at the main pul- 
sation component, complications arise when, after additional 
prewhitening s, significant components remain, or, when the 
prewhitening is ambiguous, due to the closeness of the sec- 
ondary components. There are two basic cases of comphcation: 
(i) the remnant components are not well- separated and few suc- 
cessive prewhitenings do not lead to a complete elimination of 
the secondary components; (ii) we have several well-defined 
additional peaks, whose pattern is difl^erent from the ones usu- 
ally associated with the BL phenomenon (single-sided for BLl 
and symmetrically spaced for BL2 patterns). 

Case (i) can be considered as a consequence of some non- 
stationarity in the data, leading to a nearly continuous spectral 
representation, whereas case (ii) can be related perhaps to some 
modulation of more complicated type than that of the BL phe- 
nomenon. In order to make the phenomenological classification 
of variable types as clear as possible, supported by the test to 
be presented below, we adopt the following definitions; 

PC (period changing stars); Close components (i.e., the 
ones that are close to the main pulsation frequency) can- 
not be eliminated within three prewhitenening cycles, 
or, if they can, their separation from the main pulsation 
component is less than ~ IjT, where T is the total time 
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Table 3. Definition of the signals used in the BL/PC variabihty 
test. 



Type 


Time dependence 






AM 


A(t) = A sin(nf + O); 


to(t), ip{f) 


= const. 


PM 


<p(t) = sin(n/ + Y); 


w(0,A(f) 


= const. 


FM 


aj(t) = Wo + sin(fif + T); A(t), ip{t) 


= const. 


PC 


aj(t) = wo/(l +/}t/Po); 


Ait), ipit) 


= const. 



Notes: Signal form: x{t) = (1 + A{t)) ?,m{oj{t)t + (pit)); ldq = 2n/Po; 
Po = 0.377; n = Itt/Pbl; Pel = BL period; AM=amplitude modula- 
tion; PM=phase modulation; FM=frequency modulation; PC=secular 
frequency change; O, T are arbitrary constant phases. 

span. Furthermore, except for possible harmonics, there 
are no other significant components in the spectra. 
BL (Blazhko stars): There is either only one or two sym- 
metrically spaced close components on both sides of 
the main pulsation frequency. After the second or third 
prewhitenings, except for possible harmonics, no sig- 
nificant pattern remains or the pattern is of type PC. 
MC (stars with closely spaced multiple frequency compo- 
nents): All close components are well-separated (i.e., 
frequency distances are greater than ~ l/T) and, ex- 
cept if there are only two secondary components on one 
side of the main pulsation frequency, more than three 
prewhitenings are necessary to ehminate all of them. 

We note that possible appearance of instrumental eff'ects 
(i.e., peaks at integer d"' frequencies) is disregarded in the 
above classification. Class MC is equivalent to class vM used in 
AOO and AOS. Here the new notation is merely aimed at simpli- 
fication on the occasion of the above definition. Furthermore, 
in the previous works, variables with "dominantly" BL2 struc- 
tures were included in class BL2. Now the above scheme puts 
them among the MC stars. 

Based on the following test, we note that below fre- 
quency separations less than ~ 1-5/T, distinction among the 
above types becomes more ambiguous, because in the simple 
prewhitening technique followed in this paper, the result will 
depend on the phase of the modulation. 

In order to substantiate the above definitions, here we ex- 
amine in more detail how we can distinguish between the 
PC and BL phenomena based solely on the properties of the 
prewhitened spectra. 

We use the time distribution of arbitrarily selected stars to 
generate artificial time series with modulated signal parameters 
given in Table 3. We see that the chosen types of modulation 
cover nearly all basic cases in the lowest order approximation 
(i.e., single main pulsation component with frequency ojo, lin- 
ear period change, etc.). No noise is added, because we are 
interested in differences caused by the various non-stationary 
components in the residual spectra, after prewhitened by the 
pulsation component. Nevertheless, noise plays an important 
role at low modulation levels, but a more detailed study of the 
complicated problem of non-stationary signal classification in 
the presence of noise is out of the scope of the present work. 

The same pulsation period of 0.377 d is used for 10 ran- 
domly selected stars. For each star, we scan the modulation 



frequency Q and the rate of period change yS at fixed values of 
all the other parameters. Amplitudes are adjusted to get mod- 
ulation levels that are 20% of the peak of the main pulsation 
component in the frequency spectra. Parameter /3 is changed in 
the range of (1-12) x 10"^, yielding modulation levels between 
20% and 50%. Each scan is repeated with 10 randomly selected 
phase values. By using the same code for the analysis of these 
artificial signals as the one employed for the observed data, we 
get the result shown in Fig. 6. In order to characterize the effi- 
ciency of the prewhitening, here we use the ratio A3/A1, where 
Ai and A3 denote the peak amplitudes after the first and third 
prewhitenings, respectively. Although in practice we use the 
SNR of the frequency spectra to select significant components, 
here the direct comparison of the amplitudes is more meaning- 
ful, because for noiseless signals SNR can be high even if the 
peak amplitude itself is small. 




.0004 .0008 .0012 .0016 



\Af, 



Fig. 6. Testing various modulated signals as given in Table 3. 
The ratio of the peak amplitudes after the third and first 
prewhitenings are plotted as a function of the measured fre- 
quency distance (in [d"']) from the main pulsation compo- 
nent. Dashed vertical Une shows the position of the 1/T fre- 
quency, where T denotes the total time span. Dots of different 
shades/colors are for the various modulation types: yellow/light 
gray, green/gray, red/dark gray and black are for AM, PM, FM 
and PC, respectively. The separation of simple amplitude and 
phase modulations from the more complicated ones (FM and 
PC) are clearly visible at |A/i| > 2/T. 



It is clear from the figure that PC signals remain difficult to 
prewhiten, except perhaps near and below l/T, when any types 
of signal can be prewhitened due to the proximity to the main 
pulsation component. Although periodic frequency modulation 
behaves similarly to secular period change above 1.5 /T, it can 
be prewhitened (and therefore confused) with simple ampli- 
tude and phase modulations below this limit. When the secular 
change or periodic modulation of the pulsation frequency is 
strong, they are more easily separated from simple modulated 
signals (of types AM and PM, see Table 3). This happens at 
measured modulation frequencies greater than ~ 2/T (see the 
gap starting at ~ 0.0008 d"'). Additional separation between 
AM and PM signals is observable at modulation frequencies 
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Table 6. Properties of the MDM stars 



MACHO ID 


Pi 


P2/P1 


V 


V-R 


12.10202.285 


0.398 


0.807 


18.99 


0.50 


12.10443.367 


0.337 


0.802 


18.90 


0.45 


9.4278.179 


0.327 


0.805 


18.69 


0.29 



greater than ~ 5/T , where, as expected, the AM type leaves 
the weakest trace after the third prewhitening. 

We conclude that the standard method of prewhitening we 
use for the analysis of the observed data is capable of making 
distinction between linear period change and signal modula- 
tion, but confusion occurs with periodic frequency modulations 
under observed modulation frequencies of 1.5 /J. Modulation 
types AM and FM are also difficult to distinguish, except for 
short modulation periods, when type AM can be prewhitened 
more easily. These results support our phenomenological clas- 
sification of variables given at the begiiming of this section 
even in the case of long modulation periods. 

5. New detections by using SAM 

Here we re-address the question of the incidence rate of first 
overtone BL stars by applying the method of analysis and 
scheme of classification described in the previous sections. 

The datasets used in the analysis are described in Sect. 2. 
The final statistics of classification are shown in Table 4. 
Additional details of the analysis on a star-by-star basis are 
given in Table 5. The result presented in Table 4 is based on 
the analysis of set #2 comprising the full available two-color 
data. In comparison with the similar summary of AOO (their 
Table 7), we see that the most striking difi"erence is the increase 
of the incidence rate of the BL stars by almost a factor of two. 
This change can be attributed to the following efl'ects. 

- Longer time span: +4 stars - from the comparison of the 
detections in the 'r' data. 

- The better quality and higher amplitudes in the 'b' data: 
H-28 stars - from the comparison of the set #2 'r' and 'b' 
data. 

- Extending the allowed range ofBL periods up to the length 
of the total time span: +\2 stars - based on the SAM statis- 
tics of the set #2 results. 

- Employing SAM: +9 stars - from the comparison of the set 
#2 'b' and SAM analyses. 

- False detections in AOO: -7 stars - due to the higher cutoff 
employed here for detection limit and because of the some- 
what different classification scheme. 

For the above reasons, we also find increases in the num- 
ber of other types of variables. On the other hand, three stars, 
classified earlier by AOO as RR12 have been re-classified as 
MDM, or 'mysterious double-mode'. Basic properties of these 
stars are summarized in Table 6. These objects have been re- 
classified mostly because of their somewhat extreme position in 
the color-magnitude diagram. Indeed, if we plot the derived av- 
erage magnitudes and color indices on the color-magnitude di- 
agram of Alcock et al. (2000a), the three MDM stars fall to the 



high-luminosity red edge of the region populated by RRl stars. 
We recall that in the LMC, (V) = 19.35 and {V - R) ^ 0.22, 
for the RRl stars, with V-R < 0.3 for most of the vari- 
ables (see Alcock et al. 2004). Although these objects are about 
1 mag fainter than the faintest first/second overtone double- 
mode Cepheids (see Soszyhski et al. 2000), it may still not 
be excluded that they are faint overtone Cepheids rather than 
bright and very red overtone RR Lyraes. It is also noted that 
there is a narrow overlap in the periods of these two classes of 
stars. One needs more accurate data to understand the status of 
these intriguing objects. 

A sUght change in the statistics of the above classification 
occurs if we consider double identifications due to field over- 
laps. We consider a star to be double-identified, if the following 
two conditions are satisfied simultaneously: (i) the simple dis- 
tance derived from the co-ordinates is smaller than 2 x 10"^ 
degrees; (2) the difference between the periods is smaller than 
10"^ d. We find altogether 52 double-identified objects. Among 
these there are 32 RRl-S, 8 RROl, 5 PC, 1 BLl and 3 BL2. It is 
noted that except for 3 marginal cases of PC/S ambiguities, the 
classifications of the double-identified objects are consistent. 
The multiple identifications lead to the revised incidence rates 
of 3.5% and 3.9% for the BLl and BL2 stars, respectively. We 
see that the changes are insignificant. 

5.1. Properties of the BL and MC stars 

The basic data on the 99 BL stars are summarized in Table 7. 
The description of the columns is given in the note added to the 
table. As mentioned, in the case of BL2 stars, the modulation 
frequency /bl stands for (/+ -1- /-)/2, where /+, /_ denote the 
frequencies of the A+ and A_ modulation components and /i 
the main pulsation frequency. This type of averaging is justi- 
fied, because from the 53 BL2 variables there are only 7 with 
Sf = I/+ +f_- 2/1I > 0.0001, 4 with 6f > 0.0002 and only 
one (star 15.10068.239) with 6f = 0.00048. Although this lat- 
ter value might reflect statistically significant deviation from 
equidistant frequency spacing, the majority of the the variables 
are well within the acceptable deviations due to observational 
noise (see AOO for numerical tests). 

The BL variables have been selected on the basis of the 
classification scheme described in Sect. 4. In that scheme, MC- 
type variables (i.e., those with multiple close periods) have 
been excluded from the class of BL stars. One may wonder 
if these variables could also be considered as some subclass 
of the BL variables. We think that without knowing the un- 
derlying physics behind the BL phenomenon, it is up to one's 
preference if class MC is considered as subclass BL. Although 
here we opted not to take this choice, we note that there are 
five variables among the MC stars that show symmetric peaks 
of BL2-type within the multiple peaks. In Fig. 7 we show two 
representative cases when the presence of the BL2 structure is 
obvious. 

In order to be more specific on all five stars containing 
structures of type BL2, below we give a short description of 
their frequency spectra. In the following we use the notation vq 
for the main pulsation component. 
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3.4 3.45 

Frequency 

Fig. 7. Examples of type MC variables with BL2 structures. 

Arrow indicates the position of the (prewhitened) main pulsa- 
tion component. Spectra are computed by SAM. See text for 
comments on the peak structures. 



11.8751.1740 - The frequency spectrum shows some PC-type 
remnant at vo, a symmetric pair of peaks around it and an addi- 
tional component further away (see Fig. 7). 
14.8376.548 - The frequency spectrum shows a symmetric pair 
of closely spaced peaks around vq and an additional peak fur- 
ther away (see Fig. 7). 

18.2357.757 - Classified as BL2 from the SAM analysis of 
set #1. From set #2 we found an additional close component, 
therefore the final classification has become of MC. 
19.4671.684 - The frequency spectrum shows a dominantly 
BL2 structure, but there are additional peaks spred closely. 
6.6697.1565 - The frequency spectrum shows two BL2 struc- 
tures superposed, one with large and another one with small 
frequency separations. The closer pair has also lower amph- 
tudes. 

The remaining eight stars from the MC variables show dis- 
tinct peaks without any apparent BL2-type peak structures. 
Representative examples of these variables are shown in Fig. 8. 
We see that these are indeed difi'erent from the ones generally 
classified as type BL. 

It is interesting to examine the distribution of the size and 
the position of the modulation components. For comparison, 
in Fig. 9 we show the distributions of the modulation amph- 
tudes both for the RRl and RRO BL stars. It is seen that in both 
classes there are variables with very high modulation ampli- 
tudes, exceeding 0.1 mag. This may correspond to ~ 50% and 
~ 90% relative modulation levels for the RRO and RRl stars, 
respectively. In the other extreme, at low modulation ampli- 
tudes, we find cases near 0.01-0.02 mag. From the size of the 
noise and the number of the data points, we expect about these 
amplitudes as the lowest ones to be detected in this dataset. It 
is important to note that current investigations by Jurcsik and 
co-workers (Jurcsik et al. 2005a) suggest that modulation may 
occur under 0.01 mag among fundamental mode Galactic field 
Blazhko stars. It is clear that we need more accurate data on the 
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Fig. 8. Examples of type MC variables without symmetric peak 

structures of type BL2. Arrow indicates the position of the 
(prewhitened) main pulsation component. Spectra are com- 
puted by SAM. 

Table 8. Positions of the larger modulation amplitudes 



Type 


Color 




Mo, 




BLl 


'r' 


12 


32 


0.38 




'b' 


18 


44 


0.41 




'SAM' 


21 


46 


0.46 


BL2 


'r' 


13 


22 


0.59 




'b' 


21 


46 


0.46 




'SAM' 


24 


53 


0.45 



Notes: Dataset #2 is used; N+ = number of cases when the larger 
modulation amplitude has a frequency greater than that of the main 
pulsation component; Aftot = total number of BL stars identified in the 
given color. 



LMC to check if these low modulation levels are more common 
also in other stellar systems. 

Concerning the relative positions of the modulation compo- 
nents, in Table 8 we show the number of cases obtained in the 
different datasets with the larger modulation component pre- 
ceding the main pulsation component. We see that in the case 
of RRl stars we have nearly 50% probability that this happens. 
This result is basically independent of BL type, color and ex- 
tent of the dataset. We recall that AOS derived 75% for this ratio 
for the RRO stars in the LMC. 

The present study has largely extended the range of the BL 
periods known for RRl stars. The distributions of the mod- 
ulation frequencies for the RRl and RRO stars are shown in 
Fig. 10. Although the sample is much smaller for the RRl stars, 
it is clear that there are more stars among them with short mod- 
ulation periods than among the RRO stars. Nevertheless, we 
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Fig. 9. Upper panel: distribution of the 'b' modulation ampli- 
tudes of the 99 RRl stars shown in Table 7. In the case of BL2 
stars the ampUtudes of both side peaks are included. Lower 
panel: as in the upper panel, but for the 'V modulation am- 
pUtudes of the 731 RRO stars of AOS. 



should mention that Jurcsik et al. (2005a, b) found very short 
modulation periods also among RRO stars, albeit in the Galactic 
field (i.e., SS Cnc and RR Gem with Pbl = 5.3 and 7.2 d, re- 
spectively). It is also noted that the considerable surplus in the 
RRl-BL stars with long modulation periods can be attributed 
only partially to our lower limit set on /bl- We find that among 
the 34 stars contributing to the first bin in Fig. 10 for the plot 
of RRl stars, there are only 11 with l.O/T < /bl < 1.5/7' and 
16 with 1.0/r < /bl < 2.0 /T. This implies a significant peak 
in the distribution function for /bl < 0.005 even with the omis- 
sion of these objects with very long modulation periods. The 
lack of such a peak for the RRO stars might indicate a differ- 
ence in this respect between the RRO and RRl BL stars, but, 
due to the low sample size for the RRl stars, we would caution 
against jumping to this conclusion. 

5.2. Miscellaneous variables 

In Table 9 we list the properties of the variables showing fre- 
quency spectra diflicult to classify. We see that in most cases 
the detections are barely above the noise level. Therefore, we 
cannot exclude that at least some of these variables will turn 
out to be single-periodic RRl stars, when more accurate data 
will be available. In the column "Comments" we show the pe- 
riod ratios in cases when some sort of double-mode pulsation 
can be suspected (without suggesting that those variables are 
indeed of double-mode ones). 



RR1 



q: .1 




RRO 



.05 .1 .15 



BL 



Fig. 10. Upper panel: distribution of the modulation frequen- 
cies of the 99 RRl stars shown in Table 7. Lower panel: as in 
the upper panel, but for the modulation frequencies of the 731 
RRO stars of A03. 



6. Conclusions 

This work has been motivated by an earlier result of Alcock 
et al. (2000) on the low incidence rate of the first overtone 
Blazhko RR Lyrae (RRl-BL) stars in the LMC. They derived 
a rate of 4% for these stars, which is a factor of three lower 
than that of the fundamental mode (RRO) stars (see Alcock et 
al. 2003). Soszynski et al. (2003) obtained a somewhat higher 
rate of 6% from the OGLE database. However, Alcock et al. 
(2000) used MACHO 'r' data that have somewhat lower signal- 
to-noise ratio for the objects of interest, and, in addition, the 
data analyzed then, spanned a shorter time base than the ones 
available now. This suggests that the incidence rate could be 
higher than the one deciphered earlier. In a full utilization of 
aU available data we employed a spectrum averaging method 
that enabled us to increase the detection rate by 18% in com- 
parison to the best single color (i.e., 'b') rates. 

The main conclusion of this paper is that indeed, the inci- 
dence rate of the RRl-BL stars in the LMC is higher than pre- 
viously derived from nearly the same database. Even though 
the rate is 7.5% now, it is still significantly lower than that 
of the RRO stars. This latter rate is 12%, that we also expect 
to increase slightly when the same method as the one used in 
this paper on RRl stars is to be employed also on RRO stars. 
This means, that, at least for the LMC, we can treat as a well- 
established fact that RRl-BL stars are significantly less fre- 
quent than their counterparts among the RRO stars. It is also 
remarkable that after filtering out the main pulsation compo- 
nent, the highest amplitude peak appears with equal probability 
at both sides of the frequency of the main pulsation component. 
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We recall that in the case of RRO stars there is a 75% prefer- 
ence toward the higher frequency side. The smallest modula- 
tion ampUtudes detected for RRl stars are near 0.02 mag. For 
the larger sample of RRO stars this limit goes down to 0.01 mag. 
The relative size of the modulation may reach 90% for RRl 
stars, whereas the same limit is only 50% for the RRO stars. 
Furthermore, RRl-BL stars have a relatively large population 
of short-periodic (Pel < 20 d) variables, which region is nearly 
empty for the RRO stars in the LMC. 

Since the underlying physical mechanism of the BL phe- 
nomenon is still unknown, it is also unknown if these (and 
other) observational facts will put strong (or even any) con- 
straints on future theories. Nevertheless, since these statistics 
are based on large samples, their significance is high and surely 
cannot be ignored in the theoretical discussions. For example, 
we note that the very existence of BLl stars puts in jeopardy 
the present forms of both currently available models (magnetic 
oblique rotator/pulsator by Shibahashi 2000; non-radial reso- 
nant model by Nowakowski & Dziembowski 2001; see how- 
ever Dziembowski & Mizerski 2004). 
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Table 4. Final classification of the 1354 pre-selected variables 



Classification 


Short description 


Number 


Inc. rate in RRl 


RRl-S 


Singly-periodic overtone RR Lyrae 


712 


53.5 % 


RRl-BLl 


RRl with one close component 


46 


3.5% 


RR1-BL2 


RRl with symmetric frequencies 


53 


4.0% 


RRl-MC 


RRl with more close components 


13 


1.0% 


RRl-PC 


RRl with period change 


187 


14.0 % 


RRl-D 


RRl with frequencies at integer d^' 


137 


10.3 % 


RRl-Ml 


RRl with some miscellany 


13 


1.0% 


RROl 


FU/FO double-mode RR Lyrae 


165 


12.4% 


RROl-BLl 




1 


0.1% 


RROl-PC 




5 


0.4% 


RRO-S 


Fundamental mode RR Lyrae 


1 




RRO-BLl 




1 




RR0-BL2 




1 




MDM 


Mysterious double-mode 


3 




BI 


Eclipsing binary 


16 





Table 5. Notebook of the analysis 



MACHO ID 


Period 


Type r 


Type r, b 


Type r, b 


Comments 






(AOO) 


(SAM) 


(SAM) 






[d] 


6.5 


years 


7.5 years 




10.3190.501 


0.3527946 


RRl-S 


RRl-S 


RRl-S 




10.3191.363 


0.4060709 


RROl 


RROl 


RROl 




10.3193.457 


0.2977060 


RRl-S 


RRl-S 


RRl-S 


weak peak in r at 0.4477 


10.3193.514 


0.2762716 


RRl-S 


RRl-S 


RRl-S 




10.3310.723 


0.3675948 


RRl-S 


RRl-S 


RRl-S 




10.3311.612 


0.3308571 


RRl-S 


RRl-S 


RRl-S 




10.3314.787 


0.3575062 


RRl-S 


RRl-S 


RRl-S 




10.3314.873 


0.3842321 


RRl-S 


RRl-S 


RRl-S 




10.3314.916 


0.3438987 


RRl-S 


RRl-S 


RRl-S 




10.3315.758 


0.2755685 


RRl-S 


RRl-S 


RRl-S 




10.3432.666 


0.3418437 


RRl-S 


RRl-PC 


RRl-PC 




10.3434.825 


0.2776860 


RRl-S 


RRl-S 


RRl-S 




10.3434.936 


0.2769524 


RRl-S 


RRl-S 


RRl-S 




10.3435.907 


0.3295068 


RRl-S 


RRl-S 


RRl-S 




10.3550.888 


0.3317657 


RRl-S 


RRl-S 


RRl-S 




10.3552.745 


0.2922940 


RRl-BLl 


RRl-BLl 


RRl-BLl 




10.3556.986 


0.2961494 


RRl-PC 


RRl-PC 


RRl-PC 


df=0.00038 


10.3557.1024 


0.2947668 


RRl-PC 


RRl-PC 


RR1-BL2 


df=0.00043 



Note: This table is available in its entirety at CDS ( http://cdsweb.u-strasbg.fr I. Additional data (positions, magnitudes, etc.) can be found at the 
MACHO online database I http://wwwmacho.mcmaster.ca/Data/MachoData.html^. 
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Table 7. Data on the first overtone Blazhko stars 



MACHO ID. 




Pi 




/bl 




A+ 




Ao 




A_ 


MACHO ID. 




P\ 




/bl 




A+ 




Ao 




A_ 


10.3552.745 





2922940 





037408 





0666 





1280 






6.5729.958 










074254 





1045 





1641 




" 


10.3557.1024 





2947668 





000426 





1000 





2643 





0941 


6.5730.4057 





2763214 





077598 





0662 





1249 






10.4035.1095 





3188217 





080780 





0496 





2239 






6.5850.1081 





3335893 


-0 


118341 









1617 





0326 


10.4161.1053 





2874579 





094288 





0769 





1443 






6.5971.1233 





2877179 


-0 


073476 





0380 





1868 





0666 


11.9471.780 





2859825 


-0 


113512 









2076 





0764 


6.6091.1198 





2668322 


-0 


075389 









1654 





0307 


13.5713.590 





2836353 





100165 





0636 





0870 






6.6091.877 





3206106 


-0 


038525 









1046 





0642 


13.5714.442 





3168726 





109252 





0547 





1008 






6.6094.5606 





3393854 


-0 


001244 





0593 





2476 





0637 


13.5842.2468 





2731424 


-0 


084060 









2060 





0355 


6.6326.424 





3304519 


-0 


072633 









2353 





0504 


13.5959.584 





3466458 





000993 





0286 





2629 






6.6810.616 





2883850 


-0 


000403 





0552 





2113 





0704 


13.6322.342 





2621495 





180083 





0385 





1166 






6.7054.713 





4404073 


-0 


001674 





0445 





2110 





0478 


13.6326.2765 





3304520 


-0 


072702 









2387 





0334 


6.7056.836 





4476160 





001120 





0317 





1781 





0252 


13.6810.2981 





2883863 





000419 





0691 





2050 





0659 


80.6352.1495 





3079889 


-0 


071798 









0937 





0682 


13.6810.2992 





2903975 


-0 


041503 





0909 





2153 





1318 


80.6470.2000 





3706713 





000564 





0492 





2016 





0472 


14.8495.582 





2971105 





044128 





0696 





2236 





0364 


80.6595.1599 





3323705 


-0 


036844 





0326 





1977 





0812 


14.9223.737 





3098789 





031590 





0941 





2270 





1075 


80.6597.4435 





3086445 


-0 


000950 





0368 





1503 





0519 


14.9225.776 





3551288 


-0 


054055 









2589 





0903 


80.6838.2884 





3583390 





000478 





0707 





2950 





0599 


14.9463.846 





2749961 





070745 





0470 





1293 






80.6951.2395 





3114197 


-0 


000583 





0487 





1821 





0597 


14.9702.401 





2754043 


-0 


184006 









2349 





0469 


80.6953.1751 





3517135 


-0 


028083 





0677 





2699 





0709 


15.10068.239 





2976949 


-0 


033619 





0325 





1769 





0398 


80.6957.409 





4078831 





000793 





0537 





1470 





0452 


15.10072.918 





2930774 


-0 


036838 









2971 





0649 


80.6958.1037 





2752632 


-0 


034210 





0355 





2141 





0392 


15.10311.782 





3482903 





001036 





0440 





2308 





0450 


80.7072.1545 





3583444 


-0 


001338 





0268 





1236 





0263 


15.10313.606 





2925053 





064385 





0806 





1861 





0597 


80.7072.2280 





2786360 


-0 


047201 









1486 





0359 


15.11036.255 





2869357 





078243 





1029 





1030 




" 


80.7319.1287 





2978598 





000852 





0973 





2210 





0761 


15.11280.663 





3299770 


-0 


000580 









2853 





0393 


80.7437.1678 





2781401 


-0 


088902 









1858 





1289 


18.2361.870 





3254790 





029733 





0661 





2112 





0650 


80.7440.1192 





3430250 


-0 


001029 





0582 





1660 





0663 


19.4188.1264 





3702574 





000911 





0945 





1690 





0678 


80.7441.933 





2733126 





047918 





0551 





1684 





0420 


19.4188.195 





2828193 





012127 





0683 





1752 





0665 


81.8400.901 





2739736 





006469 





0404 





1651 






2.4787.770 





3357061 


-0 


000461 





0617 





2749 





0620 


81.8518.970 





3085487 


-0 


014972 









1192 





0264 


2.5032.703 





3077973 


-0 


006788 









1864 





0344 


81.8639.1749 





3759094 





044767 





0459 





2201 





0363 


2.5148.1207 





2946848 


-0 


035316 





0611 





1824 





0664 


81.8759.832 





3038356 





074458 





0377 





2144 






2.5148.713 





3212517 


-0 


000684 





0318 





1651 





0333 


81.8879.1869 





2975649 


-0 


000735 





0614 





1871 





0660 


2.5266.3864 





2793840 





065535 





0705 





1768 




" 


82.7920.1074 





3386711 





038043 





0470 





2152 






2.5271.255 





4347510 





001169 





0658 





1821 





0541 


82.8049.746 





2987354 





071957 





0675 





2256 





0482 


2.5876.741 





2859854 


-0 


072205 









2184 





0347 


82.8286.1784 





3285773 


-0 


000494 









1893 





0357 


3.6240.450 





3836371 


-0 


001203 





0641 





2129 





0583 


82.8407.312 





3563439 





000804 





0371 





2317 






3.6597.756 





3086451 


-0 


000950 





0304 





1378 





0367 


82.8408.1002 





2962930 





035408 





0694 





2299 





0503 


3.6603.795 





2765073 


-0 


017362 









2527 





0369 


82.8525.1980 





2833982 





081435 





0355 





2378 






3.6838.1298 





3583389 





000466 





0646 





3078 





0563 


82.8526.1176 





3238862 





127004 





0321 





1293 





0312 


3.6966.427 





2972764 





078251 





0272 





2148 






82.8765.1250 





3049238 


-0 


048999 





0422 





2267 





0538 


J. /UOo.O/J 


u 




— u 


UUUd4Z 


u 




u 


Z4SJ 


u 


UO/D 


oZ.o /OD. i JUj 


u 




— u 


nn^n 1 o 
UUjUls 






u 


uyjo 


u 


UjiO 


47.1521.589 





3587236 


-0 


042303 





0575 





2303 





0634 


9.4274.644 





2636423 


-0 


122986 









1529 





0343 


5.4401.1018 





2685325 


-0 


004754 





0761 





0830 





0802 


9.4632.731 





2775501 


-0 


006687 





0310 





1220 





1233 


5.4528.1555 





2752904 


-0 


010273 









0908 





0212 


9.4636.1979 





3933065 





000434 





0254 





2201 






5.4889.1060 





2815683 


-0 


097085 









1052 





0309 


9.5117.617 





3601338 


-0 


000491 





0265 





2328 





0341 


5.5008.1902 





3774283 





073639 





0436 





1327 






9.5125.1018 





2598093 


-0 


005181 





0376 





1075 





0572 


5.5250.1501 





3551757 


-0 


000513 





0378 





2143 





0433 


9.5242.1032 





2858950 





037524 





0879 





1245 






5.5367.3432 





2854189 





029594 





0370 





1800 





0367 


9.5479.852 





3217303 





000755 





1144 





1728 





0827 


5.5368.1201 





3365130 


-0 


005552 





0287 





2049 





0554 


9.5481.746 





3581725 


-0 


000506 





0234 





2063 





0345 


5.5489.1397 





2899676 


-0 


047824 









1849 





0719 


9.5606.348 





2909644 





071233 





0226 





1168 





0209 


5.5497.3874 





2680527 


-0 


018829 









1566 





0381 

























Notes: The amplitude of the main pulsation component is denoted by Aq. The modulation amplitudes are A+ and A_, corresponding to the 
larger and smaller frequencies, respectively. Amplitudes are given in MACHO instrumental 'b' magnitudes, frequencies in [d"']. The sign of 
the modulation frequency /bl is positive, if A_ < A+ and negative, if A_ > A+. For BL2 stars /bl stands for the average of the two modulation 
frequencies (see text for further details). This table appears also at CDS jhttp://cdsvye6.u-strasbg.fr^ . 
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Table 9. List of miscellaneous stars 



MACHO ID 




Vmi 


SNR 


Source 


Comments 


10.3916.849 


3.1389921 


1.2070180 


6.61 


SAM 






11.8745.899 


3.3176002 


4.7867033 


6.79 


SAM 


vo/vmi 


= 0.6931 


1 ^ <Ti O Tl 
IJ.J / Ly. / Lj 


:>.UZzUoiU 


j.yLJjy) 1 1 


^ OA 

D.yo 




Vo/VmI 


= 0.7718 


18.2717.787 


2.8643249 


4.9297051 


6.55 


SAM 


vo/vmi 


= 0.5810 


18.3202.956 


2.2855398 


5.1908177 


6.96 


r 






19.3575.541 


3.3698536 


0.5856270 


7.91 


SAM 






3.7451.484 


3.3048194 


2.7630866 


7.85 


b 


vmi/vo 


= 0.8361 


80.6348.2470 


2.9571188 


1.2069597 


6.65 


SAM 






80.7558.650 


2.7828759 


0.8240633 


6.99 


r 






81.8519.1395 


2.9463890 


0.8372214 


6.96 


r 






9.5004.750 


3.2877325 


0.1604569 


7.42 


SAM 






9.5122.363 


3.1033541 


4.1679579 


6.82 


SAM 


Vo/Vm/ 


= 0.7446 


9.5241.382 


3.1777911 


1.3826455 


6.79 


SAM 







Notes: Vq denotes the frequency (in [d ']) of the main pulsation component, Vmi stands for the frequency of the miscellaneous component with 
SNR given in the next column. Interesting period ratios are given in the conmient lines. 



